To bind or not to bind: The //-dibaryon in light of chiral effective field theory 
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I. INTRODUCTION AND SUMMARY 
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We analyse the quark mass dependence of the binding energy of the ff-dibaryon in the framework 
of chiral effective field theory. We show that the SU(3) breaking effects induced by the differences of 
the pertinent two-baryon thresholds (AA, EAT, EE) have a very pronounced impact that need to be 
incorporated properly in future lattice QCD simulations. We also point out that if the i/-dibaryon 
is a two-baryon bound state, its dominant component is EA rather than AA, which is a consequence 
I ■ of the approximate SU(3) flavor symmetry of the two-baryon interactions. 

1^) PACS numbers: 13.75.Ev,12.39.Fe;14.20.Pt 
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• In 1977 Jaffe predicted a deeply bound 6-quark state with / = 0, J = 0, S = — 2 from the bag-model, called 
| the .ff-dibaryon |l|. Subsequently, many experimental searches for the .ff-dibaryon were carried out, but so far no 
convincing signal was found 0]. However, recently evidence for a bound ff-dibaryon was claimed based on lattice 
QCD calculations Extrapolations of the simulations, performed for pion masses > 400 MeV, to the physical 

mass suggest that the .ff-dibaryon could be either loosely bound or move into the continuum 0, H| . 

In this paper, we analyze various issues related to these lattice results in the framework of chiral effective field 
Q-i, theory (EFT) for the baryon-baryon {BB) interaction at leading order (LO) in the Weinberg counting. This scheme 
has proven successful for the few data on hyperon-nucleon scattering Q and also for the bounds that exist on the 
BB interactions with strangeness S = —2 10]. There is one low-energy constant (LEC) in the S = —2 sector, 
corresponding to the SU(3) flavor-singlet channel, that can not be fitted by hyperon-nucleon data and can be fine- 
tuned to produce a bound H with a given binding energy. This framework also allows us to study the quark mass 
dependence 1 of the binding energy of the ff , see the calculations of the quark mass dependence of the deuteron binding 
energy in (TlT[l4l | . Another important issue to be addressed here is how this quark mass dependence is affected when 
the SU(3) breaking manifested in the masses of the octet baryons is accounted for. In the real world the BB thresholds 
are not degenerate but they all differ for the relevant AA, EE and EN systems and that has very definite dynamical 
consequences, as we will show in what follows. 

The pertinent results of our study can be summarized as follows: 



(i) We have analyzed the effective range expansion for the A A So channel, assuming a loosely bound H dibaryon. 
It shows a very different behaviour to the case of the deuteron in the neutron-proton 3 Si channel. In fact, any 
attraction supplied by the flavor-singlet channel contributes with a much larger weight to the HiV interaction 
than to AA, according to SU(3) flavor symmetry, so that the .ff-dibaryon should predominantly be a SAT bound 
state. 

(ii) We observe that, for pion masses below 400 MeV, the dependence of the binding energy of the ff is linearly 
decreasing with decreasing pion mass, in agreement with the findings in Ref. [7J. In particular, a H binding 
energy (BE) adjusted to the value found by NPLQCD 0] at M„ = 389 MeV, is reduced by 7 MeV at the physical 
pion mass. For larger pion masses, this dependence is weakened. Note, however, that for such large pion masses 
this should only be considered a trend as the chiral EFT is constructed for masses/momenta well below the 
chiral symmetry breaking scale. 

(iii) We find a much more drastic effect caused by the SU(3) breaking related to the values of the three thresholds 
AA, EE and SAT. For physical values the BE of the H is reduced by as much as 60 MeV as compared to a 
calculation based on degenerate (i.e. SU(3) symmetric) BB thresholds. Translating this observation to the 
situation in the HAL QCD calculation, we see that the bound state has disappeared at the physical point. 
For the case of the NPLQCD calculation, a resonance in the AA system might survive. 



1 Because of the Gell-Mann-Oakes-Renner relation, the pion mass squared is proportional to the average light quark mass. Therefore, the 
notions "quark mass dependence" and "pion mass dependence" can be used synonymously. 



2 



5 = 



Channel 



AW -> NN 



AN -> A AT 
A AT -> EAT 
SAT ->■ EAT 
EAT ->■ EAT 



AA 
AA - 
AA - 
EN - 
EN 

EE - 
EN - 
EN 
EA - 
EA - 
EE - 



AA 
HAT 
EE 

■ EN 

■ EE 

EE 

■ EN 

■ EA 
EA 
EE 

■ EE 



Isospin 



C 2 



27 +c ,S 3) 



-C 27 + C 8s ) 



ij {C 27 + 9C 8 



10 
C 27 



^ (27C* 27 + 8C 8s +5C 1 ) 

(18C* 27 - 8C* 8s - 10C 1 ) 
^§ (-3C 27 + 8C 8s - 5C 1 ) 
-L (12C 27 + 8C* 8s + 20 C 1 ) 
*f (2C 27 + 8C 8s - 10C* 1 ) 
i (C 27 + 24C 8s +15C 1 ) 
| (2C 27 + 3C 8 =) 
^(C" 27 -C 8 =) 
i (3C 27 + 2C 8 *) 

C 27 



Isospin 







C, 



si 



C 1 



2 

Q10 



C 8 " +c 
-c 8 » + c 



10* 
10 



•) 




TABLE I: Various LO baryon-baryon contact potentials for the 1 So and 3 Si partial waves in the isospin basis. C 27 etc. refers 
to the corresponding SU(3)f irreducible representation. 



Our manuscript is organized as follows: In Sec. UU we recall the basic formalism of the BB interaction in the 
framework of chiral EFT. Sec. Mil contains a detailed discussion of the quark mass dependence of the BE of the H 
and the influence of the SU(3) breaking through the various two-baryon thresholds. In Sec. IIVI we try to make direct 
contact to the results published by the NPLQCD and HAL QCD collaborations. 



II. THE BARYON-BARYON INTERACTION TO LEADING ORDER 



For details on the derivation of the chiral BB potentials for the strangeness sector at LO using the Weinberg power 
counting, we refer the reader to Refs. 0,113 » see also Refs. Here, we just briefly summarize the pertinent 

ingredients of the chiral EFT for BB interactions. 

The LO potential consists of four-baryon contact terms without derivatives and of one-pseudoscalar-meson ex- 
changes. The LO SU(3) f invariant contact terms for the octet BB interactions that are Hermitian and invariant 
under Lorentz transformations follow from the Lagrangians 

C 1 = C} (B a B b (T t B) b (TiB) a ) , C 2 = Cf (B a (T t B) a B b (T t B) b ) , 

£ 3 = Cf (B a (TiB) a ) (B b (TiB) b ) . (1) 

Here a and b denote the Dirac indices of the particles, B is the irreducible octet (matrix) representation of SU(3) f , 
and the Tj are the usual elements of the Clifford algebra [9] . As described in Ref. [9J , to LO the Lagrangians in Eq. ([T]) 
give rise to only six independent low-energy coefficients (LECs), the Cfs in Eq. JT]), due to SU(3) f constraints. They 
need to be determined by a fit to experimental data. It is convenient to re-express the BB potentials in terms of the 
SU(3)f irreducible representations, see e.g. Refs. [l8|,[l9j]. Then the contact interaction is given by 

V=^{l-a x - <j 2 ) C 1S0 + 1(3 + CTl ■ <r 2 ) C 3S1 , (2) 

and the constraints imposed by the assumed SU(3) f symmetry on the interactions in the various BB channels for the 
1 5'o and 3 Si partial waves can be readily read off from Table |U 

The lowest order SU(3) f invariant pseudoscalar-meson-baryon interaction Lagrangian embodying the appropriate 
symmetries was also discussed in [9(. The invariance under SU(3) f transformations implies specific relations between 
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the various coupling constants, namely 
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Here / = gA/^F^, where is the nucleon axial- vector strength and is the weak pion decay constant. We use the 
values g A = 1.26 and F w = 92.4 MeV. For a, the F/(F + £>)-ratio [1 . we adopt the SU(6) value: a = 0.4, which is 
consistent with recent determinations of the axial- vector coupling constants [20]. 

The spin-space part of the LO one-pseudoscalar- meson-exchange potential is similar to the static one-pion-exchange 
potential in chiral EFT for nucleon-nucleon interactions, see e.g. [2 11 ] (recoil and relativistic corrections give higher 
order contributions), 

v b 1 b 2 ^b[b' 2 f f (gi ' q) ( (T 2 • q) 

V - -Ib^pIb^p q2 + M 2 ' W 

where Mp is the mass of the exchanged pseudoscalar meson. The transferred and average momentum, q and k, are 
defined in terms of the final and initial center-of-mass (cm.) momenta of the baryons, p' and p, as q = p' p and 
k = (p' + p)/2. In the calculation we use the physical masses of the exchanged pseudoscalar mesons. The explicit 
SU(3) breaking reflected in the mass splitting between the pseudoscalar mesons and, in particular, the small mass of 
the pion relative to the other members of the octet leads to sizeable differences in the range of the interactions in the 
different channels and, thus, induces an essential dynamical breaking of SU(3) symmetry in the BB interactions. The 
r] meson was identified with the octet r\ (rjg) and its physical mass was used. 

The reaction amplitudes are obtained from the solution of a coupled-channels Lippmann-Schwinger (LS) equation 
for the interaction potentials: 



The label v indicates the particle channels and the label p the partial wave. is the pertinent reduced mass. The 



on-shell momentum in the intermediate state, q v , is defined by y/s = Uvn^ + q% + y?^ 2 ^ + 9 2 .- Relativistic 

kinematics is used for relating the laboratory energy Tj a b of the hyperons to the cm. momentum. 

In [l(| the LS equation was solved in the particle basis, in order to incorporate the correct physical thresholds. 
Since here we are only interested in the H dibaryon we work in the isospin basis. Then for J — 0, I = 0, and S = —2 
we have to consider the three coupled channels AA, EN and EE. We use the following (isospin averaged) masses 
to a = 1115.6 MeV, m s = 1192.5 MeV, m H = 1318.1 MeV, and m N = 939.6 MeV so that the AA, EN, and EE 
thresholds are at 2231.2, 2257.7, and 2385.0 MeV, respectively. Furthermore, the potentials in the LS equation arc 
cut off with a regulator function, exp [— (p' 4 + p 4 ) /A 4 ], in order to remove high-energy components of the baryon 
and pseudoscalar meson fields I22| . We consider cut-off values in the range 550, 700 MeV, similar to what was 
used for chiral NN potentials |22| . 

The imposed SU(3) flavor symmetry implies that only five of the six LECs contribute to the YN interaction, 
namely C 27 , C 10 , C 10 , C 8s , and C 8a , cf. Table Q] These five contact terms were determined in [9| by a fit to the YN 
scattering data. Since the NN data cannot be described with a LO EFT, SU(3) constraints from the NN interaction 
were not implemented explicitly. As shown in Ref. @ , a good description of the 35 low-energy YN scattering can 
be obtained for cutoff values A = 550, 700 MeV and for natural values of the LECs. The sixth LEC (C 1 ) is only 
present in the S = —2 channels with isospin zero, cf. Table [I] There is scarce experimental information on these 
channels that could be used to fix this LEC, but it turned out that the quality of the existing data do not really allow 
to constrain its value reliably [l0| . Even with the value of the sixth LEC chosen so that C^g^ AA — 0, agreement 
with those data can be achieved. In this case a scattering length of a^ n — —1.52 fm [l(| is obtained. Analyses of the 
measured binding energy of the double-strange hypernucleus A j{He [23| suggest that the AA scattering length could 
be in the range of -1.3 to -0.7 fm [3S]. 



III. QUARK MASS DEPENDENCE OF THE BINDING ENERGY AND SU(3) BREAKING EFFECTS 



Chiral effective field theory itself does not allow one to make any predictions with regard to the existence of the 
H dibaryon because, as stated, one of the contact terms (C 1 ) occurs just in the channel in question, cf. Table |H 
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and has to be determined from there. If we take over the values for C 27 and C 8s as fixed from the YN data and 
assume that C 1 = then we find no bound state for AA - neither in LO [H nor based on the preliminary NLO results 
[2T| . However, if we assume that C 1 7^ and vary its value then in both cases a near-threshold bound state can be 
produced for values of natural size. 

At the same time, the framework of chiral effective field theory in which our AA interaction is derived is very well 
suited to shed light on the general characteristics of a ii-dibaryon, should it indeed exist. In particular, it allows us 
study the implications of the imposed (approximate) SU(3) flavor symmetry and to explore the dependence of the 
properties of an assumed ii-dibaryon on the masses of the relevant mesons and baryons. The latter aspect is a rather 
crucial issue in view of the fact that the available lattice QCD calculations were not performed at the physical masses 
of the involved particles. 

To start our discussion, let us assume that the ii-dibaryon is a (loosely) bound BB state and that its binding energy 
Eh is similar to that of the deuteron D. For pedagogical purposes we fix the value of the flavor-singlet LEC C 1 in 
such a way that = jd = 0.23161 fm {E = — 7 2 /tob, where tub is either tojv or toa), because of the well-known 
relation between the binding energy and the effective range parameters [H, [29[ 

1 1 2 

a 2 

This relation is very well fulfilled for the deuteron and the corresponding neutron-proton 3 Si scattering length (a — 
5.43 fm) and effective range (r — 1.76 fm). One would naively expect that the same should happen for the H- 
dibaryon. However, it turns out that the corresponding results for AA in the So partial wave are quite different, 
namely a = 3.00 fm and r = —4.98 fm. Specifically, the effective range is much larger and, moreover, negative. Clearly, 
the properties of the ii-dibaryon are not comparable to those of the deuteron, despite the fact that both bound states 
are close to the elastic threshold. Indeed, if one recalls the expressions for the relevant potentials as given in Table U 

^AA^AA = 1 , 27C 27 + 8C 8 S + ^ ySN^SN = 1 ,^27 + 8C <8 S + 

40 v ' 40 v ; 

one can see that the attraction supplied by the SU(3) flavor-singlet state (C 1 ) contributes with a much larger weight 
to the EN channel than to AA. This indicates that the presumed ii-dibaryon could be predominantly a EN bound 
state. We have confirmed this conjecture by evaluating explicitly the phase shifts in the AA and EN channels, cf. the 
discussion in the next section. Indeed, one finds that the phase shift for the EN channel is rather similar to the NN 
3 Si case. Specifically, the EN ( 1 S'o) phase shift 5(qsN) fulfills <5(0) — <5(oo) = 180° in agreement with the Levinson 
theorem. The AA ( 1 5'o) phase behaves rather differently and satisfies <5((/aa = 0) — S(oo) = 0. Note that there have 
been earlier discussions on this issue in the context of S = —2 baryon-baryon interactions derived within the quark 
model 0, IM|. 

Let us now consider variations of the masses of the involved particles. The dependence of the H binding energy on 
the pion mass is displayed in Fig. Q] (left). For the case considered above, enlarging the pion mass to around 400 
MeV (i.e. a value corresponding to the NPLQCD calculation Q) increases the binding energy to around 8 MeV and 
a further change of M n to 700 MeV (corresponding roughly to the HAL QCD calculation [![) yields then 13 MeV, cf. 
the solid line. Readjusting C 1 so that we predict a H binding energy of 13.2 MeV for M n = 389 MeV, corresponding 
to the latest result published by NPLQCD [f|, yields the dashed curve. Fig. Q] includes also results of a calculation 
where C 1 was fixed in order to reproduce their earlier value of around 16 MeV Q (dash-dotted curve), to facilitate 
a more direct comparison with chiral extrapolations [7|, [8j based on that value. It is obvious that the dependence on 
Mjr we get agrees well - at least on a qualitative level - with that presented in Ref . Q . Specifically, our calculation 
exhibits the same trend (a decrease of the binding energy with decreasing pion mass) and our binding energy of 9 
MeV at the physical pion mass is within the error bars of the results given in 0- On the other hand, we clearly 
observe a non- linear dependence of the binding energy on the pion mass. As a consequence, scaling our results to the 
binding energy reported by the HAL QCD collaboration @] (30-40 MeV for M n m 700 - 1000 MeV) yields binding 
energies of more than 20 MeV at the physical point, which is certainly outside of the range suggested in Ref. [?}■ 
However, we note that for such large pion masses the LO chiral EFT can not be trusted quantitatively. We remark 
that in our simulations the curves corresponding to different binding energies remained roughly parallel even up to 
such large values as suggested by the HAL QCD collaboration. 

Our finding that any ii-dibaryon is very likely a bound EN state rather than a AA state, which follows from the 
assumed (approximate) SU(3) symmetry of the interaction, suggests that not only the pion mass but also the masses 
of the baryons play a significant role for the concrete value of binding energy. Indeed, the physical difference between 
the AA and EN thresholds of around 26 MeV implies that the ii-dibaryon considered above is, in reality, bound by 
roughly 28 MeV with respect to its "proper" threshold. Accordingly, one intuitively expects that in a fully SU(3) 
symmetric case, where the masses of all octet baryons coincide, the bound state would remain basically fixed to the 
EN threshold and then would lie also about 28 MeV below the AA threshold. In the concrete case of J = 0, I = 0, 
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FIG. 1: Dependence of the binding energy of the if-dibaryon on the pion mass Mn (left) and on the E mass ms (right). The 
solid curve correspond to the case where C 1 is fixed such that Eh = — 1.87 Me V for physical masses while for the dashed curve 
C 1 is fixed to yield E H = -13.2 MeV for AU = 389 MeV. The pion mass dependence is also shown for E H = -16.6 MeV 
(dash-dotted curve). The asterisks and crosses represent results where, besides the variation of ros, m~ + mjv = 2mA is 
assumed so that the HiV threshold coincides with that of the AA channel. The vertical (dotted) lines indicate the physical AA 
and EE thresholds. 



S = — 2 we are dealing with three coupled channels, namely AA, HiV, and EE. Since we know from our experience 
with coupled-channel problems [1, I32T - I341 ] that coupling effects are sizeable and the actual separation of the various 
thresholds plays a crucial role we investigated also the dependence of the H binding energy on the thresholds (i.e. on 
the E, and on the H and N masses). Corresponding results are displayed in Fig. [T]on the right side. The solid curve 
is again the result based on our reference case with a binding energy of -1.87 MeV for physical masses of the pion 
and the baryons. When we now decrease the E mass so that its nominal threshold of 2385 MeV moves downwards 
and finally coincides with the one of the AA channel (2231.2 MeV), we observe a rather drastic change in the H 
binding energy. Note that the direct interaction in the EE channel is actually repulsive for the low-energy coefficients 
C 27 and C 8s fixed from the YN data plus the pseudoscalar meson exchange contributions with coupling constants 
determined from the SU(3) relations Eq. And it remains repulsive even for C 1 values that produce a bound 
ii-dibaryon. But the coupling between the channels generates a sizeable effective attraction which increases when 
the channel thresholds come closer. The dashed curve is a calculation with the contact term C 1 fixed to simulate the 
binding energy (-13.2 MeV) of the NPLQCD collaboration at M n = 389 MeV. As one can see, the dependence on the 
binding energy on the E mass is rather similar. The curve is simply shifted downwards by around 4.5 MeV, i.e. by the 
difference in the binding energy observed already at the physical masses. The asterisks and crosses represent results 
where, besides the variation of the EE threshold, the EN threshold is shifted to coincide with that of the AA channel. 
This produces an additional increase of the H binding energy by 20 MeV at the physical EE threshold and by 9 MeV 
for that case where all three BB threshold coincide. Altogether there is an increase in the binding energy of roughly 
60 MeV when going from the physical point to the case of baryons with identical masses. This is significantly larger 
than the variations due to the pion mass considered before. Note that we have kept the pion mass at its physical 
value while varying the BB thresholds. 



IV. APPLICATION TO LATTICE QCD RESULTS 



After these exemplary studies let us now try to connect with the published H binding energies from the lattice 
QCD calculations [3, 4J. The results obtained by the HAL QCD collaboration are obviously for the SU(3) symmetric 
case and the corresponding masses are given in Table I of Ref. Q. Thus, we can take those masses and then fix the 
LEC C 1 so that we reproduce their H binding energy with those masses. To be concrete, we use m ps — 673 MeV and 
m,B = 1485 MeV, and fix C 1 so that Eh = —35 MeV. When we now let the masses of the baryons and mesons go to 
their physical values the bound state moves upwards, crosses the AA threshold, crosses also the S7V threshold and 
then disappears. In fact, qualitatively this outcome can be already read off from the curves in Fig. [T]by combining 
the effects from the variations in the pion and the baryon masses. Based on those results one would expect a shift of 
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FIG. 2: Phase shifts for AA ( 1 So) (left) and SA^ ( So) (right) as a function of the pertinent laboratory energies. The solid line 
is the result for our reference BB interaction that produces a bound H at Eh = —1.87 MeV. The other curves are results for 
interactions that are fine-tuned to the H binding energies found in the lattice QCD calculations of the HAL QCD (dashed) 
and NPLQCD (dash-dotted) collaborations, respectively, for the pertinent meson (pion) and baryon masses as described in the 
text. 



the H binding energy in the order of 60 to 70 MeV for the mass parameters of the HAL QCD calculation. 

In case of the NPLQCD calculation we take the values provided in Ref. [Hj|, i.e. rriN = 1151.3 MeV, m\ = 
1241.9 MeV, m s = 1280.3 MeV, and m H = 1349.6 MeV. Those yield then 17 MeV for the EN-AA threshold separation 
(to be compared with the physical value of roughly 26 MeV) and 77 MeV for the EE-AA separation (physical value 
around 154 MeV). With those baryon masses we fix again the LEC C 1 so that we reproduce the H binding energy 
given by the NPLQCD collaboration, namely E H = -13.2 MeV @. We use also M n = 389 MeV, but we take the 
physical masses for the other pseudoscalar mesons (K and 77). Again we let the masses of the baryons and of the pion 
go to their physical values. Also here the bound state moves upwards and crosses the AA threshold. However, in the 
NPLQCD case the state survives and remains below the EN threshold at the physical point. Specifically, we observe 
a resonance at a kinetic energy of 21 MeV in the AA system or, more precisely, a quasi-bound state in the EN system 
around 5 MeV below its threshold. 

The results reported so far were all obtained with the LECs as fixed from the YN data for the cutoff value 
A = 550 MeV [9 1 ]- In order to investigate the stability of our results we considered also variations of the cutoff 
mass A in the LS equation Eq. (J5J) between 550 and 700 MeV, as in and repeated the exercise described above. 
Those variations led to changes of the predicted resonance energy in the order of 1 MeV only. Even an exploratory 
calculation that utilizes the low-energy coefficients determined in a preliminary NLO study of the YN system [27| 
yielded practically the same result. 

It is interesting to observe that the chiral extrapolation of the lattice QCD results performed by Shanahan et al. 
Q yields results that are similar to ours. In that reference the authors conclude that the ii-dibaryon is likely to be 
unbound by 13±14 MeV at the physical point. Let us emphasize, however, that our values are not really comparable 
with theirs. In our analysis we assume that the iJ-dibaryon is actually a bound BB state - which seems to be the 
case also in the lattice QCD studies 0, H|. On the other hand, in Ref. [§J it is assumed that the H is a compact, 
multi-quark state rather than a loosely bound molecular state. How such a genuine multi-quark state would be 
influenced by variations of the BB thresholds is completely unclear. It depends, among other things, on whether and 
how strongly this state couples to the AA, EN, and EE channels. So far there is no information on this issue from 
lattice QCD calculations. Clearly, in case of a strong and predominant coupling to the AA alone, variations of the 
EE and EN would not influence the H binding energy significantly. However, should it couple primarily to the EN 
and/or EE channels then we expect a sensitivity to their thresholds comparable to what we found in our study for 
the case of a bound state. 

Finally, for illustrative purposes, let us show phase shifts for the 1 So partial wave of the AA and EN channels. 
This is done in Fig. [2j The solid line is the result for our reference BB interaction that produces a loosely bound 
H dibaryon with Ejj = —1.87 MeV. The phase shift for the EN channel (left side) is rather similar to the one for 
the 3 5i NN partial wave where the deuteron resides, see e.g. (22[. Specifically, it starts at 180°, decreases smoothly 
and eventually approaches zero for large energies, fulfilling the Levinson theorem. The result for AA ( 1 S'o) (left side) 
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behaves rather differently. This phase commences at zero degrees, is first negative but becomes positive within 20 
MeV and finally turns to zero again for large energies. The dashed curve corresponds to the interaction that was 
fitted to the result of the HAL QCD collaboration and reproduces their bound H dibaryon with their meson and 
baryon masses. The phase shift of the SiV channel, calculated with physical masses, shows no trace of a bound state 
anymore. Still the phase shift rises up to around 60° near threshold, a behavior quite similar to that of the 1 Sq NN 
partial wave where there is a virtual state (also called antibound state [3(1]). Indeed, such a virtual state seems to be 
present too in the S7V channel as a remnant of the original bound state. The effect of this virtual state can be seen in 
the AA phase shift where it leads to an impressive cusp at the opening of the SiV channel, cf. the dashed line on the 
left side. In the EN phase shifts for the NPLQCD case (dash-dotted curve) the presence of a bound state is clearly 
visible. The corresponding AA phase shift exhibits a resonance-like behavior at the energy where the (quasi-bound) 
H dibaryon is located. 

Given the present uncertainties of the lattice QCD results for the H dibaryon binding energy we have focused here 
primarily on qualitative features of the H that can be inferred from chiral effective field theory. Clearly, once more 
precise lattice data become available one should also perform a careful assessment of the uncertainties involved in the 
EFT calculation. 
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